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Section  1 


INTRODUCTION  AND  SUMMARY 


This  report  documents  the  results  and  conclusions  of  the  Add-On  to 
Contract  No.  F33615-76-C-1031.  The  goal  in  this  effort  was  to  show  the 
feasibility  of  constructing  a nonlinear  lens  capable  of  imaging  in  the 
infrared  wavelength  spectrum.  Except  for  wavelength  of  operation,  this 
lens  has  optical  characteristics  similar  to  the  existing  visual  spectrum 
nonlinear  lens  which  was  developed  during  a previous  contract1*. 

That  lens  takes  advantage  of  the  "variable  acuity"  characteristic  of 
human  vision  to  reduce  the  amount  of  information  that  must  be  transmitted  in 
a wide  f ield-of-view  high  resolution  imaging  system.  A brief  description  of 
this  concept  is  presented  in  Appendix  A.  A laboratory  brassboard  of  that 
system  was  constructed  under  contract  to  the  Office  of  Naval  Research  (ONR)2 
and  is  presently  under  final  testing  at  MCAIR. 

The  system  developed  for  ONR  operates  in  the  visual  wavelength  spectrum 
and  consequently  is  usable  only  during  daylight  hours.  The  goal  of  this 
effort  was  to  determine  the  feasibility  of  achieving  the  same  wide  field-of- 
view,  high  resolution  capability  during  night-time  hours  by  operating  the 
sensing  lens  in  the  infrared  wavelength  spectrum.  If  successful,  this  sensor 
system  could  fully  support  human  vision  under  clear  night  conditions;  a capa- 
bility impossible  to  achieve  with  conventional  FLIR  imaging  systems. 

A brief  description  of  this  effort  and  the  results  are  now  presented.  A 
flow  diagram  of  the  overall  effort  is  shown  in  Figure  1.  It  consists  of  two 
phases,  an  analysis  phase  and  a design  phase.  In  the  analysis  phase,  the 
optimum  parameters  and  constraints  are  determined  that  are  necessary  for  the 
actual  lens  design  effort  which  will  be  accomplished  in  the  design  phase. 

A computer  controlled  optimization  analysis  indicated  that  the  lens 
should  be  designed  as  an  F/2.0  with  an  optical  quality  represented  by  a point 
spread  function  in  the  focal  plane  which  has  a one  sigma  width  (radius)  of 
16  microns.  The  best  operating  region  was  found  to  be  8-14  microns. 

A narcissus  analysis  shown  on  Figure  1 was  performed  on  the  optical  relay 
early  in  the  effort  to  assist  in  establishing  element  curvatures  and  spacings 
of  optical  elements  close  to  the  focal  plane.  As  the  actual  lens  design 
evolved,  the  narcissus  analysis  was  expanded  to  include  those  elements. 

The  lens  design  phase  proved  conclusively  that  an  infrared  nonlinear 
lens  is  feasible.  It  can  be  built  with  the  required  optical  characteristics 
(F/2.0,  16  micron  blur,  8-14  micron  wavelength  spectrum)  and  is  of  reasonable 
size  (8  inch  clear  aperture  maximum) . Narcissus  images  were  found  to  be  well 
below  the  detector  Noise  Equivalent  Power  (NEP)  when  current  generation  scan- 
ners are  used.  However,  when  used  with  second  generation  monolithic  focal 
plane  array  detectors,  dynamic  electronic  compensation  of  detector  outputs  to 
reduce  narcissus  images  to  an  acceptable  level  will  be  required  because  of 
the  larger  cold  assembly  area  and  anticipated  better  sensitivity. 

* Number  shown  as  superscript  refers  to  the  corresponding  item  in  the  List  of 
Refe rences. 
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FIGURE  1 

ANALYSIS  FLOW  DIAGRAM 
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Section  2 


ANALYSIS 


The  analysis  portion  of  this  effort  proceeded  along  the  lines  shown  on 
Figure  1.  Basically  it  is  divided  into  two  separate  efforts,  the  Optimization 
Analysis  and  the  Narcissus  Analysis.  They  were  run  in  parallel  which  was 
dictated  by  the  tight  schedule  for  the  supply  of  data  required  by  the  lens 
design  effort. 

The  key  portion  of  the  analytic  effort  is  the  optimization  effort.  As 
shown  on  the  figure  this  required  the  definition  of  all  inputs,  the  develop- 
ment of  computer  programs,  the  performance  of  size/performance  parameter 
tradeoffs,  and  integration  of  results  to  facilitate  selection  of  the  best  non- 
linear lens  design  parameters. 

The  narcissus  analysis  paralleled  the  optimization  analysis  so  that  guide- 
lines for  minimizing  narcissus  could  be  established  for  the  lens  design  effort 
concurrent  with  the  availability  of  other  optical  design  parameters.  Again 
refering  to  the  flow  diagram  of  Figure  1,  this  meant  that  the  narcissus  theory 
had  to  be  developed,  since  none  was  available  in  sufficient  detail.  The  goal 
of  the  first  narcissus  effort  was  to  assess  the  contribution  of  the  optical 
relay  since  this  impacted  on  the  rear  element  design  of  the  nonlinear  lens. 
After  the  lens  design  began  to  evolve,  narcissus  analysis  was  run  on  the  lens 
to  determine  its  contribution  to  narcissus  image  generation.  Only  minor 
changes  to  lens  design  to  reduce  narcissus  was  possible  because  of  time  and 
effort  constraints.  Details  of  these  analyses  are  presented  below. 

2.1  OPTIMIZATION  ANALYSIS 

APPROACH 

The  overall  objective  of  this  effort  is  to  establish  the  parameters  for 
an  infrared  foveal  lens  design.  Since  specific  mission  requirements  and/or 
state-of-the-art  in  infrared  variable  acuity  technology  are  not  defined  at 
this  time,  the  approach  selected  was  a computerized  optimization  which  allowed 
interactive  variation  of  key  parameters  while  size  and  performance  were 
assessed.  An  optimum  size/performance  balance  was  achieved  and  a lens  was 
preliminarily  designed  for  these  optimum  parameters.  A flow  diagram  of  this 
approach  is  shown  in  Figure  2. 

Optimization  philosophy  is  as  follows:  Experience  from  our  visual  lens 
work  shewed  that  lens  size  (linear  dimensions)  increase  inversely  with  both 
optical  quality  (blur)  and  F/number.  Conversely  performance  measured  by  Mini- 
mum Resolvable  Temperature  (MRT)  improves  directly  with  these  two  parameters. 
The  approach  shewn  on  Figure  2 was  developed  to  determine  the  best  performance 
but  with  reasonable  size. 

Inputs  for  the  study  were 

(a)  Parameters  from  our  existing  visual  spectrum  lens 
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(b)  Present  and  future  scanner  and  detector  parameters* 

(c)  System  inputs 
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FIGURE  2 

OPTIMIZATION  ANALYSIS  TECHNICAL  APPROACH 


In  general,  these  inputs  fall  within  three  categories.  These  are  constant 
inputs  that  do  not  change  in  the  study,  variable  inputs  which  were  changed 
discretely  at  various  points  during  the  studies,  and  tradeoff  parameters 
that  were  tested  over  a range  of  values  to  determine  the  effect  of  each  on 
performance  and  size.  The  latter  are  F/number  and  blur  (previously  discussed) 
and  the  number  of  effective  scan  lines  per  picture  height  (N) . This  addi- 
tional tradeoff  parameter  surfaces  early  in  the  study  as  a key  parameter 
because  of  its  major  influence  on  system  resolution  and  the  possibility  of  its 
control  in  future  detector  development.  The  variable  inputs  are  spectral 
wavelength  regions  and  scanner/detector  parameters.  The  spectral  wavelength 
bands  used  in  the  study  are  in  the  3-5  micron  region  and  8-14  micron  regions. 
Since  these  regions  require  different  detectors,  D*  and  electronic  detector 


response  were  changed  when  spectral  band  was  changed.  The  other  variable 
inputs  are  the  scanner  parameters  which  are  established  by  present  and  future 
| technology.  A state-of-the-art  serial  scanner  and  a parallel  scanner  were 

evaluated  separately3.  H^ever,  the  primary  emphasis  in  this  study  was 
placed  on  the  utilization  of  a second  generation  monolithic  focal  plane  array 
scanner.  Such  a scanner  was  theorized  in  both  spectral  bands  while  the  state- 
1 of-the-art  scanners  were  studied  in  the  long  wavelength  band  only.  The  func- 

tion of  these  inputs  and  the  rationale'  for  their  selection  in  the  analysis 
are  described  below. 

LENS  OPTICAL  PARAMETERS 

Starting  at  the  upper  left  of  Figure  2,  the  optical  requirements  of  dis- 
tortion, field-of-view,  and  fabrication  constraints  are  used  to  define  focal 
length  and  image  size  requirements  of  the  nonlinear  lens.  The  basis  for  the 
selection  of  these  inputs  is  as  follows:  (1)  The  amount  of  distortion  used 
was  made  equal  to  the  existing  lens.  This  function  provides  a good  match  to 
human  visual  acuity  and  is  believed  to  be  near  the  maximum  distortion  tech- 
nically possible.  In  addition,  any  sensing  system  constructed  with  this 
function  will  be  compatible  with  the  display  systems  constructed  on  related 
contracts.  (2)  A field-of-view  (FOV)  goal  of  160°  was  utilized  for  the  study. 
This  value  was  selected  because  no  significant  reduction  in  either  size  or 
bandwidth  is  achieved  for  smaller  90°-100°  fields-of-view.  In  addition  this 
FOV  does  not  impose  any  difficulty  in  design.  (3)  Fabrication  of  our  existing 
visual  spectrum  lens  has  shown  that  the  rear  element  spherics  are  extremely 
difficult  to  produce  if  the  diameter  is  less  than  0.75  inch. 

The  relationships  established  by  these  parameters  are  shown  in  Figure  3. 
These  curves  were  used  to  select  an  on-axis  focal  length  of  2 inches.  This 
then  established  an  image  size  of  0.72  inch  for  160°  coverage.  A curve  for  a 
100°  FOV  is  also  shown  on  Figure  3 to  emphasize  the  point  made  above  that 
little  is  gained  through  FOV  reduction. 


Image  Height  - in. 
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FIGURE  3 

FOCAL  LENGTH/IMAGE  SIZE  PARAMETRICS 

*A11  analyses  and  results  pertaining  to  current  technology  scanners  was 
extracted  from  a related  NAVY  contract.3 
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OPTICAL  RELAY  PARAMETERS 

When  the  image  size  is  known,  it  is  possible  to  define  the  optical  relay 
that  couples  the  image  to  the  detector  plane.  The  basic  geometry  for  this  is 
shown  in  Figure  4.  The  optical  relay  must  be  defined  because  it  establishes 
how  the  nonlinear  lens  F/number  appears  to  the  detector  (Figure  2)  for  sensi- 
tivity calculations  and  to  establish  how  the  detector  dimensions  translate  to 
the  lens  focal  plane  (Figure  3),  required  for  spatial  frequency  calculations. 
The  required  ratio  (Lc/Ld)  is  the  ratio  of  the  nonlinear  lens  image  size  (H) 
to  the  (effective  detector  scan  plane  height  (Hjj)  (Figure  4).  The  lens  image 
size  is  C . 36  inch.  The  effective  detector  height  Hd  depends  on  the  detector 
assembly  and  scan  technique  used. 
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FIGURE  4 

OPTICAL  RELAY  PARAMETERS 


In  general  this  height  is: 

N a 

Hn  = n 

D 2 ovsc 

where 

N = the  effective  number  of  scan  lines 
a = vertical  detector  dimension 

y 

novsc  = line  overscan  ratio 


(1) 
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Since  the  optics  are  collimated  between  the  two  lenses 


'c  LD 


Substituting  Equation  (2)  and  (1)  and  rewriting 


N a n 


The  detector  size  (Ah  ) measured  in  the  y direction  and  in  the  lens  image 
plane  is  F 

Ah  a 

_Z  = JL  (4) 

L Ln 

c D 

Using  the  focal  length  ratio  of  Equation  (3)  and  substituting  into  (4)  results 

in  . 

SL  q u 

y c 2 H 


Note  that  Ah  is  independent  of  detector  size.  In  the  x direction,  the 
detector  siz£  is 

L o „ a 

Ah  - a -£  - ■■  2 H -* 
x x Ln  N rj  a 
D ovsc  y 

If  a and  a are  related  by  a constant  (K) 
x y J 


un  %T  — rv  an  ' ' ' 

x N n y 

ovsc 

Since  these  increments  are  not  dependent  on  detector  size,  the  detector  spatial 
ffTF’s  will  also  be  independent  of  detector  size. 

The  F/number  of  the  detector  (F#D)  in  terms  of  the  nonlinear  lens  (F//) 
will  new  be  developed  from  Figure  4.  The  ray  bundle  diameter  (D)  shown  on 
Figure  4 is  defined  by 

L 

D - ~ (8) 

The  F/nunfcer  (F^)  of  the  nonlinear  lens  is  also  equal  to  the  F/number  of  the 
relay  (F^c) . The  F/number  of  the  detector  (F|(D)  is 
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(9) 


Solving  for  D and  setting  Equation  (9)  equal  to  (8)  results  in 


(10) 


The  F/number  ratio  in  terms  of  scanner  parameters  may  be  defined  by  substitut- 
ing Equation  (3)  into  (10) 


Nan 
y ovsc 

2 H 


(ID 


An  interesting  observation  can  be  derived  from  this  equation.  From  Reference 
(4)  signal  power  falling  on  a single  detector  element  is  directly  proportional 
to  detector  area  (Ap)  and  inversely  to  detector  F/number  (F#D) , viz. 


and  since  the  detector  area  is 


(12) 


a a 


x y 


(13) 


Therefore 


P 


D 


(14) 


If  Equation  (11)  is 
F/number  (F^)  which 


P 


D 


used  to  express  signal  power  in  terms  of  nonlinear  lens 
is  by  definition  also  relay  lens  F/number  (F„  ) 

''  /"i 


4 H 


ovsc 


(15) 


If  the  iiiiage  height  (H)  and  overscan  ratio  (n)  are  constant  we  may  write  a 
proportionality 


P 


D 


2 2 
N F/ 


(16) 
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Therefore  signal  pcwer  is  not  dependent  on  detector  area  but  only  on  F/number 
and  scan  line  count.  In  NVL3  the  noise  power  is  proportional  to  the  square 
| root  of  detector  area.  The  result  is  that  the  NEDT  and,  therefore,  MRT  are 

directly  proportional  to  the  square  root  of  detector  area.  Therefore,  the 
MRT  is  directly  proportional  to  detector  linear  dimension  if  the  detector 
remains  symmetrical  as  its  size  varies.  If  the  number  of  scan  lines  and 
nonlinear  lens  F/number  remain  constant,  smaller  detectors  result  in  better 
performance.  Since  a simple  linear  relationship  exists  between  performance 
and  detector  size,  it  was  not  retained  as  a tradeoff  parameter  for  second 
generation  detector  analysis.  Instead,  the  smallest  detector  dimension  that 
was  practical  was  used,  viz. 


a = a = 0.001  inch 
x y 

This  leaves  F//  and  N as  the  only  scanner  related  tradeoff  parameters.  F // 
was  varied  from  F/1.0  to  F/11.0  and  scan  lines  from  131  to  2000.  For  current 
generation  scanners,  detector  sizes  and  N are  fixed.  Thus  only  F //  is  a trade- 
off parameter.  These  fixed  values  of  detector  size  and  N are  shown  in 
Figure  5. 
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FIGURE  5 

CURRENT  GENERATION  SCANNER  PARAMETERS 
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New  the  overscan  ratio  will  be  defined.  The  overscan  ratio  (h^g,,)  is  the 
ratio  of  detector  vertical  subtense  (Ah  ) to  actual  scan  line  spacing  and  may 
be  defined  by  the  equation:  ^ 


a L N 
y c 

^ovsc  2H 


(17) 


To  put  a value  on  this  parameter,  the  type  of  scanning  to  be  used  must  be 
theorized.  For  the  primary  tradeoff  analysis  a second  generation  hybrid  scan 
technique  was  used,  i.e.,  a parallel  array  with  serial  time  delay  and  integra- 
tion (TDI) . For  these  systems  an  overscan  ratio  of  one  was  used.  For  the 
current  generation  scanners,  the  serial  scanner  used  a ratio  of  1.5  while  the 
parallel  scanner  utilized  a ratio  equal  to  one. 
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T 


In  the  second  generation  analyses,  the  number  of  serial  detectors  was 
always  equal  to  nine.  This  number  was  based  on  customer  inputs  which  stated 
that  the  payoff  of  using  more  than  nine  was  questionable  when  such  things  as 
switching  noise  and  relative  image  detector  motion  during  TDI  were  considered. 
The  current  generation  serial  scanner  also  employed  nine  TDI  detectors. 

The  remaining  scanner  parameters  of  D*  was  selected  as  follows:  For  3-5 
micron  systems,  InSb  as  a detector  was  assumed  with  a 

D*  = 5 x 1011 


For  the  8-14  micron  system  the  same  material  was  selected.  The  D*  value  is 


D*  = 5 x 10 
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These  values  include  improvements  due  to  effective  cold  shielding.  Manufac- 
turer's actual  improvement  data  was  used  to  estimate  the  improvement  rather 
than  theoretical  values.  Detector  equivalent  electronic  response  BW  was 
determined  from  the  manufacturer's  literature.  This  parameter  was  found  to  be 
significant  only  in  the  case  of  the  serial  scanner. 

LENS  SIZE  PARAMETRICS 

During  initial  phases  of  the  study,  the  lens  size  parametrics,  i.e,,  size 
as  a function  of  F/number  and  optical  blur  of  the  existing  visual  spectrum 
lens  was  used.  These  parametrics  are  shown  in  Figure  6.  As  actual  lens 
design  progressed,  these  parametrics  were  updated.  Hwever,  this  was  not 
found  to  be  necessary.  All  optimized  designs  fit  these  data  very  accurately. 


PERFORMANCE  CALCULATIONS 

System  performance  in  terms  of  MRT  was  computed  using  the  Night  Vision 
Laboratory  (NVL)  FLIR  model  which  was  suitably  modified  for  the  VARVS  system. 
This  model  and  its  modifications  is  described  in  Appendix  B.  Other  than 
parameters  previously  discussed,  this  model  required  the  parameters  listed  on 
Figure  7.  The  basis  for  selection  of  the  more  salient  parameters  is  listed 
below. 

(a)  Wavelength  Range  - For  both  spectrums  these  were  selected  as  regions 
of  good  atmospheric  transmission  over  reasonable  path  lengths. 

(b)  Atmospheric  transmission  - Original  plans  were  to  use  0.5  but  were 
changed  to  1.0  to  reflect  system  performance  rather  than  the 
extremely  variable  mission  performance. 

(c)  Optical  transmission  - Based  on  state-of-the-art  coating  techniques, 
a 0.5%  reflectance  for  each  of  10  elements  (two  surfaces  each). 
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FIGURE  6 


SIZE  PARAMETRICS  FOR  A 2-IN.  FOCAL  LENGTH  LENS 

(d)  (e)  System  electronic  response  - Since  the  NVL  model  includes  no  internal 
noise  sources  other  than  the  detector,  electrical  response  was  made 
insignificant  for  all  systems  by  assuming  flat  response  out  to  the 
highest  spactial  frequency  (fx)  °f  interest  which  was  determined  by 

f > — , where  I is  the  instantaneous  FOV.  This  assumption  is 

X 21 FOV  F0V 

entirely  adequate  for  the  NVL  model  because  (1)  response  is  limited 
by  detector  geometry,  and  (2)  noise  bandwidth  is  limited  by  human 
observer  factors  such  as  eye  integration  time  and  the  matched  filter 
function  (see  Appendix  B) • 

(f)  Display  blur  - The  characteristics  of  the  direct  view  display  were 
used  for  this  parameter.  This  unit  has  a blur  sigma  in  angular  terms 
of  0.26  minutes  arc.  This  is  about  10  times  better  than  the  antici- 
pated system  resolution.  When  utilizing  the  NVL  model  it  was  most 
convenient  to  convert  this  to  linear  spatial  blur  at  the  same  focal 
length  as  the  nonlinear  lens  which  is  3.9  microns. 

(g)  Display  magnification  - The  VARVS  concept  inherently  has  unity 
magnification. 


(h)  Frame  rate  - A TV  compatible  frame  rate  of  30  frames /second  is  uti- 
lized for  all  analyses  with  two  fields  per  frame  (2/1  interlace). 

(i)  Eye  integration  time  - The  NVL  recommended  value  of  0.2  sec.  was 
used  for  this  parameter. 

(j)  Temperature  - The  usual  300°K  value  was  used  for  nominal  scene 
temperature. 
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(k)  Emissivity  - An  emissivity  of  0.8  is  utilized  for  both  target  and 
background. 

(l)  Threshold  S/N  - The  NVL  value  of  2.25  was  used  for  this  parameter. 

(m)  Vibration  sigma  - The  vibration  sigma  selected  represents  the 
residual  sensor  motion  after  stabilization.  A value  of  1 micron  was 
used  and  this  does  not  cause  any  significant  performance  loss. 

PARAMETRIC  ANALYSIS  AND  SELECTION  RATIONALE' 

Using  the  parameters  defined  above,  the  analysis  and  selection  proceeds 
as  follows.  For  each  wavelength  band  and  scanner  type  the  tradeoff  parameters 
are  iterated  throughout  their  ranges.  For  each  set  of  tradeoff  parameters, 

MRT  vs  fx  and  lens  size  are  printed.  Selection  then  proceeds  as  follows. 

First  a specific  value  of  fx  is  chosen  which  is  the  spatial  frequency  at  which 
the  optimization  will  be  made.  Then  for  each  scan  line  value  (N) , the  data 
is  screened  to  determine  the  remaining  tradeoff  parameters  (F/number  and  blur) 
that  yield  the  smallest  system  for  a range  of  MRT  performance  levels.  This 
procedure  is  shown  in  Figure  8.  The  spatial  frequency  is  then  varied  and  the 
entire  process  is  repeated.  This  procedure  was  selected  because  clear  size 
minimums  were  observed  as  a function  of  N,  similar  to  those  shown  on  Figure  8. 
These  minimums  were  recorded  along  with  their  system  parameters  and  used  to 
plot  size  vs  MRT  for  each  spatial  frequency  of  optimization  as  shown  on  the 
right  of  Figure  8.  These  data  were  used  to  locate  the  point  where  a reduction 
of  MRT  is  not  warranted  by  the  corresponding  increase  in  size.  The  knee 
points  of  the  curves  represent  the  best  selection.  Beyond  these  points, 
excessive  size  results  with  little  improvement  in  MRT.  Finally  the  MRT 
curves  for  the  selected  systems  (one  optimum  for  each  design  fx)  are  plotted 
as  shown  on  the  lower  right  of  Figure  8.  A subjective  choice  is  then  made  as 
to  which  is  most  acceptable  from  a mission  performance  and  vehicle  installa- 
tion standpoint.  The  results  of  the  analysis  will  now  be  presented. 

RESULTS  OF  OPTIMIZATION  ANALYSIS  FOR  SECOND  GENERATION  DETECTOR  TECHNOLOGY 

The  second  generation  scanner  size  parametrics  as  a function  of  effective 
scan  lines  are  shown  in  Figures  9 through  12  for  the  3-5  micron  spectrum  and 
in  Figures  13  through  16  for  the  8-14  micron  spectrum.  Spatial  frequencies 
of  optimization  are  0.3,  0.5,  0.7  and  1.0  cycle /mr.  Separate  figures  are  pre- 
sented for  each  of  these  frequencies.  Performance  MRT  values  range  from 
0.05°C  to  0.5°C.  These  ranges  of  both  MRT  and  spatial  frequency  were  chosen 
because  they  cover  the  anticipated  requirements.  Some  comments  on  these  data 
are: 

(a)  The  3-5  micron  curves.  Figures  9 through  12  are  presented  in  order 
of  increasing  spatial  frequency  of  optimization.  Note  how  the 
optimum  scan  minimum  size  increases  with  spatial  frequency  of 
optimization  with  clear  minimums  occurring  on  Figures  11  and  12 
(0.7  and  1.0  cycle/mr.) 

(b)  At  any  particular  spatial  frequency  of  optimization,  minimums 
occur  at  the  same  scan  frequency  regardless  of  the  MRT  requirement. 


Program 
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Nominal  Wavelength  Band 
Actual  Wavelength  Band  Used 
Detector  Type 
D* 

Detector  Width 

Detector  Height 

Half-Image  Height  on  Nonlinear 
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Nonlinear  Lens  On-Axis  Focal  Length 
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Focal  Length  Ratio 

Relay  Focal  Length 


Relay  F/Number 

F/Number  Evaluated 

Optics  Transmission 
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Detector  Time  Response  BW 

Electronics  BW 

Atmospheric  Transmission 
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Eye  Integration  Time 
Total  Scan  Lines  Evaluated 
Vertical  Scan  Efficiency 
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Horizontal  Scan  Efficiency 
Overscan  Ratio 
Frame  Rate 
Display  Spot  Size  (la) 


DSTAR 
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microns  Sigma 


FSTAR 


2nd  Gen  Det  System 


State  of  Art  Systems 


microns  SIGMAD 


XN  is  the  effective  number  of  scan  lines,  9 is  the  relay  field-of-view  angle 


FIGURE  7 


PARAMETERS  USED  IN  SYSTEM  EVALUATION 
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• Compute  MRT  (fx,  Fno,  a,  N) 

/F#n°n  Nn  Size] 

I fx  MRT  I 


°1  N1  Size 
P#2  a1  ^1  Size 

°1  ^1  Sizel 


fx  MRT 


• For  Selected  fx's  Plot 
Size  vs  N for  Specific 
MRT  values* 


fx3 

fx2 


fXl 


N 


• Plot  Minimum  Size 

System  for  Each  fX/MRT 
Combination.  Select  Best 
Lens  Parameters  Based  on 
ASize/AMRT 


•This  step  not  necessary  for  current  scanners  where  N is  specified. 


Plot  MRT  Curves 
for  Selected 
Optimum  Systems 


FIGURE  8 

PARAMETRIC  ANALYSIS  AND  SELECTION 
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FIGURE  9 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  3-5  n Frequency:  0.5  c/mr 
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FIGURE  10 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  3-5  m Frequency:  0.3  c/mr 
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FIGURE  11 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  3-5 Frequency:  0.7  c/mr 
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FIGURE  12 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  3-5  m Frequency:  I.Oc/mr 

(c)  Trends  in  the  8-14  micron  data  on  Figures  13  through  16  are  approxi- 
mately the  same  as  for  the  3-5  micron  case.  However,  sizes  are 
significantly  lower  for  the  same  spatial  frequency.  This  is  clearly 
illustrated  at  the  next  step  in  optimization  shown  on  Figure  17. 

In  Figure  17,  line  rate  is  eliminated  as  a variable  by  selecting  the  systems 
from  the  previous  figures  that  have  the  optimum  line  rate  (that  which  yields 
minimum  size  for  the  same  performance).  This  curve  is  used  to  make  the  final 
tradeoffs  between  spectral  regions  and  between  size  and  performance.  The 
former  is  clearly  seen  by  comparing  performance  of  systems  of  similar  sizes. 

For  example,  if  0.7  c/mr  spatial  resolution  is  desired,  an  8 inch  diameter 
system  would  have  a MRT  of  0.12°C  in  the  8-14  micron  spectral  region  and  MRT 
of  0.3°C  in  the  3-5  micron  spectral  region.  Another  way  of  making  the  compari- 
son is  on  a constant  MRT  performance  basis.  For  example,  if  0.1°C  MRT  is 
required  at  0.7  cycle/mr,  a 3-5  micron  system  would  have  a diameter  of  almost 
14  inches  while  the  8-14  micron  would  have  8.5  inch  diameter.  Similar  trends 
exist  at  all  performance  and  size  levels.  Therefore,  the  8-14  micron  spectral 
region  seems  to  be  a clear  choice. 

The  curves  of  Figure  17  can  also  be  used  to  locate  the  best  size  or  per- 
formance combination.  This  occurs  at  the  knee  of  the  curves  where  size  starts 
increasing  excessively  while  performance  improves  only  slightly.  These  points 
are  shown  on  Figure  18  for  the  six  optimum  systems  (three  spatial  frequencies, 
two  spectral  regions).  A description  of  these  optimum  systems  is  listed  in 
Figure  19.  A direct  comparison  of  the  performance  of  these  six  systems  is  pre- 
sented in  Figure  20.  Here  the  MRT  vs . spatial  frequency  is  shown  for  each  sys- 
tem. These  curves  show  again  the  size  advantages  of  the  8-14  micron  region  and 
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FIGURE  13 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  8-14  pt  Frequency:  0.3  c/mr 


Effective  TV  Lines 

FIGURE  14 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  8-14 /i  Frequency:  0.5  c/mr 
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FIGURE  15 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  8-14  n Frequency:  0.7  c/mr 


FIGURE  16 

2nd  GENERATION  DETECTOR  LENS  SIZE  REQUIREMENTS 

Wavelength:  8-14  n Frequency:  I.Oc/mr 
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FIGURE  17 

SIZE/PERFORMANCE  PARAMETRICS  2ND  GENERATION  DETECTOR  TECHNOLOGY 

3-5  n and  8-14  p Spectral  Regions 


FIGURE  18 

OPTIMUM  SYSTEMS  OF  THE  2nd  GENERATION  DETECTOR  TECHNOLOGY 

3-5 1 u and  8-14  n Spectral  Regions 


19 


Optimization  Parameter 

Performance 

Size 

System  Parameters 

Spectral 

Spatial 

MRT 

(°C) 

Clear 

Effective 

Blur 

Region 

Frequency 

Aperture 

Scan 

F/No. 

Diameter 

(microns) 

(cycles/mr) 

Dia  (in.) 

Lines 

(microns) 

8-14 

■ I 

4.0 

365 

2.8 

22 

8-14 

HI 

481 

2.0 

16 

8-14 

1.0 

pm 

650 

1.4 

12 

0.140 

5.6 

365 

2.0 

22 

WM 

9.8 

481 

1.4 

18 

1.0 

20.5 

650 

1.0 

12 

GP77  0292  18 


FIGURE  19 

PARAMETERS  OF  OPTIMUM  SYSTEMS 

show  how  each  system  performs  at  spatial  frequencies  other  than  that  of  design 
optimization.  The  latter  point  is  the  most  significant.  Note  that  the  system 
optimized  for  0.7  cycles/mr  actually  is  capable  of  performing  well  beyond 
1.0  cycle/mr.  This  brings  up  an  interesting  point.  How  does  performance  of 
an  optimum  system  vary  with  scan  line  count  in  regions  above  the  frequency  of 
optimization?  Scan  line  count  could  conceivably  affect  performance  adversely 
since  it  was  selected  on  the  basis  of  a relatively  low  spatial  frequency  com- 
pared to  the  ultimate  system  performance.  To  study  this,  the  8-14  micron 
system  optimized  for  0.7  cycles/mr  was  evaluated  at  all  scan  line  rates. 
Results  are  shown  on  Figure  21.  Note  that  the  selected  system  (481  lines) 
yields  the  best  performance  of  any  line  rate  in  the  range  of  0.4  to  1.0  cycle/ 
mr.  Beyond  this  spatial  frequency,  650  lines  gives  slightly  better  per- 
formance. Beyond  1.13  cycles/mr,  the  937  line  system  produces  slightly 
better  performance.  All  other  line  rates  degrade  performance  significantly. 

It  is  concluded  that  the  optimized  system  yields  nearly  optimum  performance 
at  all  spatial  frequencies. 

Performance  drop  at  lew  scan  line  count  is  obviously  because  the  detector 
spatial  MTF  becomes  dominant.  The  cause  of  the  performance  drop  at  high  scan 
line  count  is  not  so  obvious.  It  is  caused  by  a loss  in  detector  sensitivity 
because  of  the  smaller  detector  area.  This  can  be  seen  most  easily  by  observ- 
ing the  image  of  the  detector  array  at  the  nonlinear  lens  image  plane.  Since 
this  image  is  always  the  same  size,  more  scan  lines  must  be  closer  spaced,  and 
hence  smaller  in  area. 

Final  selection  is  based  on  size  and  performance  requirements.  For  the 
high  performance  manned  aircraft  application  at  which  this  system  is  directed, 
a size  in  the  range  of  8-10  inches  is  acceptable.  This  size  is  dictated  by 
window  swept  volume  considerations  and  aerodynamic  considerations.  Referring 
to  Figure  17  or  18,  this  dictates  systems  optimized  for  0.7  cycle/mr  are 
required.  Either  3-5  micron  or  8-14  micron  systems  fit  the  size  requirement 
and  can  easily  support  human  visual  performance.  However,  the  8-14  micron 
system  is  significantly  smaller.  For  these  reasons,  the  8-14  micron  system 
is  believed  to  be  the  best  choice.  Therefore,  lens  design  for  the  second 
generation  detector  technology  should  be  directed  towards  achieving  an  F/2.0 
lens  with  a maximum  blur  sigma  of  16  microns. 
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FIGURE  20 

MRT  CURVES  FOR  OPTIMUM  SYSTEMS 

2nd  Generation  Detector  Technology 
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FIGURE  21 

MRT  FOR  OPTIMUM  0.7  c/mr  SYSTEM  WITH  2ND  GENERATION  DETECTOR 
AS  A FUNCTION  OF  EFFECTIVE  SCANNING  LINE  COUNT 
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RESULTS  OF  OPTIMIZATION  ANALYSIS  FOR  CURRENT  GENERATION  SCANNER  TECHNOLOGY 


Current  generation  scanner  technology  analysis  was  performed  under  ONR 
contract1.  Since  the  line  rate  is  defined  for  the  current  generation  scanners, 
selection  proceeds  immediately  to  Size  vs  MRT  curves,  similar  to  Figure  17  and 
which  is  presented  on  Figure  22.  Note  that  sizes  are  similar  to  those  of  the 
second  generation  scanner  systems  but  the  performance  (MRT)  is  not  as  good  as 
the  second  generation  systems.  The  F/number  and  blur  of  the  optimum  systems 
for  0.7  cycle/mr  are  approximately  the  same  as  they  were  in  the  second  genera- 
tion system.  These  are  noted  on  Figure  22.  Based  on  this  observation,  a non- 
linear lens  having  an  F/2.0  and  a blur  of  16  microns  maximum  was  also  selected 
as  the  best  lens  for  current  generation  scanners. 

ANALYSIS  CONCLUSIONS 

The  performance  of  all  systems  are  compared  on  Figure  23.  The  advantages 
of  the  second  generation  detector  technology  is  readily  apparent  from  this 
data.  The  performance  of  the  current  generation  scanners  are  comparable  but 
the  serial  scanner  has  somewhat  better  performance. 

As  a final  step,  system  MTF's  are  plotted  for  all  four  scanner  systems. 
These  are  presented  separately  in  Figures  24  through  27.  Close  observation 
of  the  MTF  curves  shews  the  relative  balance  achieved  between  all  significant 
degrading  phenomenon. 

The  analysis  indicates  that  the  F/2.0  lens  with  an  optical  blur  sigma  of 
16  microns  operating  in  the  8-14  micron  spectral  region  is  optimum.  This  lens 
can  be  used  effectively  with  and  is  near  optimum  for  either  current  or  second 
generation  detectors. 

2.2  NARCISSUS  ANALYSIS 

Two  types  of  narcissus  analysis  were  conducted  during  this  effort.  The 
first  was  conducted  to  devise  the  most  efficient  method  of  relaying  the  lens 
image  to  the  detector  focal  plane.  In  the  second  effort,  actual  narcissus 
imagery  were  located  and  its  effect  on  system  operation  was  predicted.  Details 
of  these  analyses  are  presented  below.  Basic  narcissus  theory  is  presented 
in  Appendix  C. 

RELAY  NARCISSUS 

Concern  with  the  possible  requirement  for  a field  lens  at  the  nonlinear 
lens  focal  plane  necessitated  an  early  parametric  analysis  of  the  optical 
relay  with  respect  to  narcissus  effects.  Narcissus  is  caused  by  the  cold 
detector  "seeing"  itself  as  a result  of  reflections  from  relay  optical  sur- 
faces. This  relay  is  shown  in  general  form  in  Figure  28.  The  function  of  this 
relay  is  to  reimage  the  lens  image  of  the  real  world  onto  the  detector.  This 
would  be  an  easy  optical  task  if  it  were  not  for  the  scanning  mechanism  which 
essentially  translates  the  real  world  image  across  the  detector.  The  scanning 
mechanism  requires  that  all  ray  bundles  cross  at  a common  plane  shown  as  the 
Effective  Scanner  stop  of  diameter  D on  Figure  28.  A problem  arises  because 
the  optical  output  of  the  nonlinear  lens  is  not  compatible,  with  this  require- 
ment. The  chief  rays  emerging  from  the  image  plane  are  not  "well  behaved"  as 
they  are  in  a conventional  lens  and  therefore  will  not  arrive  at  a common 
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FIGURE  23 

DETECTOR  PERFORMANCE  AS  A FUNCTION  OF  RESOLUTION 
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Response 


SERIAL  DETECTOR/ NONLINEAR  SYSTEM  MTFs 


F/2.0  Blur  Circle  16m  488  Line  8- 14p  Spectral  Region 


FIGURE  25 

PARALLEL  DETECTOR/FOVEAL  LENS  SYSTEM  MTF’S 

F/2.0  Blur  Circle  16p  480  Lines  8-14  Spectral  Region 
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Response 


Detector  Time  Response, 
Electronics,  and 
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FIGURE  26 

2nd  GENERATION  OETECTOR/NONLINEAR  LENS  SYSTEM  MTF's 
F/1. 4 Blur  Circle  18/i  481  Line  3-5p  Spectral  Region 


FIGURE  27 

2nd  GENERATION  DETECTOR/NONLINEAR  LENS  SYSTEM  MTF'S 
F/2.0  Blur  Circle  16p  481  Line  8- 14p  Spectral  Region 
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a)  Center  Scan 


OPrr-OM-i 

FIGURE  28 

RVS  RELAY  AND  NARCISSUS  GEOMETRY 

plane  after  collimation.  This  necessitates  use  of  an  aspheric  field  lens  at 
the  image  plane.  While  field  lenses  create  no  special  problem  in  visual  sys- 
tems, they  can  cause  very  serious  "Narcissus"  problems  in  the  infrared  spec- 
trum. The  geometry  of  narcissus  is  also  shown  in  Figure  28. 

It  is  easy  to  see  how  optical  surfaces  near  any  image  plane  generate  the 
worst  possible  narcissus  effects  because  they  reflect  and  return  the  image  back 
through  the  system  with  maximum  efficiency.  The  field  lens  is  in  an  ideal 
location  to  create  this  problem.  For  this  reason,  the  narcissus  theo  y 
developed  in  Appendix  C was  applied  to  this  relay.  For  this  study,  a reflec- 
ing  surface  was  theorized  with  curvature  (R)  and  separation  (S  ) from  the 
nominal  image  location. 

Narcissus  theory  requires  that  the  location,  size  and  strength  of  the 
narcissus  image  be  known.  This  requires  definition  of  the  distance  (S)  from 
the  narcissus  image  to  the  image  plane,  magnification  of  the  image  (M) , 
and  reflectance  of  the  reflecting  surface  (t),  as  shown  in  Figure  28.  The 
first  two  can  be  determined  from  geometric  optical  theory.  Using  the  defini- 
tion of  distances  shewn  on  Figure  28,  the  distance  (S)  and  magnification  (M) 
are: 


S = 


c 


LD  (S2-S1> 

+ <VS1>  (L_Lc"Ld) 


(18) 
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where 


S = — 

2 2S1  + R 

and  the  magnification  is: 


M = - 


Lc  S2 

S1  [Lt2  + (Sj-Sj)  <L-Lo-Ld)  J 


The  narcissus  image  power  on  an  on-axis  detector  of  area  (A,)  is: 

d 

x y 

Pn  = — y in  watts 

U S 


where 


•t  = Surface  reflectance  determined  from  optical  material  and 
coating  characteristics 

f)Q  pay 

x = (-xz — or  M — r— ) whichever  is  smaller 
2LD  2 

DS  aCay 

y = (-rr-  or  M — r*~)  whichever  is  smaller 
zld  i 

a = Effective  cold  assembly  width,  height 

C3x ) cay 


This  is  the  worst  case  value  which  will  occur  at  the  center  display 
scan  position  as  shown  in  Figure  28.  For  other  scan  positions,  narcissus 
returns  are  much  more  difficult  to  compute.  However  the  maximum  scan  angle 
beyond  which  the  narcissus  return  vanishes  is  relatively  easy  to  predict. 
This  occurs  at  the  angle  (6  ) shewn  in  Figure  28  where  the  marginal  ray  from 
the  detector  strikes  the  curved  reflecting  surface  normally.  At  this  scan 
angle  only  this  ray  returns.  At  greater  scan  angles  the  reflection  does  not 
pass  through  the  stop  and  thus  no  narcissus  exists.  This  maximum  scan  angle 
is 


@L  = arc  tan 


u 

r ^ i 

+ Lc  - A. 


or  image  position 


Sff-R  + Lc  * A] 
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If  the  image  position  (x)  is  divided  by  the  format  width,  the  portion  of  the 
display  occupied  by  the  narcissus  image  can  be  established. 
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The  above  theory  was  applied  to  the  four  scanners  postulated  in  this 
study.  Characteristics  of  these  scanners  and  their  required  relay  components 
are  listed  in  Figure  29.  Relay  focal  lengths  were  computed  to  match  the 
detector  plane  to  the  lens  focal  plane  as  described  previously  in  Section  2.1. 
The  scanner  input  geometry  was  scaled  to  accommodate  the  F/2.0  lens  output  and 
to  provide  a full  1 inch  by  1 inch  scan  format.  The  detector  parameters  of 
these  four  scanners  that  are  required  for  narcisuus  calculations  are  also 
listed  in  Figure  29.  The  data  for  the  first  three  scanners,  the  Worst  Case 
Serial,  Probable  Serial,  and  Parallel,  is  taken  from  ONR3. 


Scanner  Type 

Effective  Cold 
Assembly 
Dimensions  cm 

7 

NEP  - watts 

Lens  Focal 
Length  - in. 

Separation 

(in.) 

(L) 

1 

acax 

acay 

Aq  • cm* 

— 

Collimat- 

ing 

<LC> 

Detector 

(Lq) 

1 

Worst  Case  Serial 

0.3513  0.234 

12.5x10“® 

8.05x10“ 11 

1.291 

1 

4 

1 

Probable  Serial 

0.116 

0.0066 

12.5x10“® 

8.05x10“ 11 

1.29 

1 

4 

1 

Parallel 

0.0087  1.882 

2.5x10“® 

2.57x10“n 

2.46 

2.67 

4 

2 

2nd  Generation 

0.038 

1.504 

12.5x10— 6 

3.8x10“12 

1.29 

1 

4 

1 

GP77-032S-28 


FIGURE  29 

SCANNER  PARAMETERS 


The  modified  serial  scanner  includes  a reflective  detector  surround  to 
redvce  narcissus.  This  narcissus  reduction  technique  is  also  employed  in  the 
parallel  scanner  and  is  theorized  for  all  second  generation  scanners. 

The  equations  developed  in  previous  paragraphs  were  applied  to  the  relay 
designed  for  each  scanner  type  to  determine  how  the  strength  and  the  size  of 
the  narcissus  return  varies  with  respect  to  the  image  location  (S^)  as  loca- 
tion of  the  field  lens  surface  and  its  radii  of  curvature  (R^)  are  changed. 
Separation  (S^)  was  varied  + 2 cm  and  four  radii  were  used,  (+  1 cm,  + ») . 

The  results  are  presented  on  Figures  30  to  37.  Figures  30  and  31  show  nar- 
cissus for  the  worst  case  serial  scanner  with  positive  and  negative  field 
lens  radii.  Note  on  Figure  30  that  with  a flat  field  lenses  (infinite  curvature 
Rl  = °°) , separation  must  be  maintained  greater  than  approximately  + 0.5  cm  to 
make  narcissus  acceptable  (below  the  NEP).  If  the  field  lens  has  positive 
radii,  the  negative  separation  must  be  increased.  As  shown  on  Figure  31,  nega- 
tive radii  on  the  field  lens  reverses  this  situation. 
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Separation,  Si  - cm 


FIGURE  30 

SERIAL  SCAN  ASSEMBLY 

Positive  Case 
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Power  - Watts  x 10  Raster  Ratio 


a^*!^  0.3513  cm  r ^ * — 1 cm 

acay  = 0.234  cm  — — — • R,  = -“cm 


-2-10  1 2 

Separation,  S-|  - cm  gp77-«m-3 

FIGURE  31 

SERIAL  SCAN  ASSEMBLY 

Negative  Case 


31 


— t*-* 


FIGURE  32 

SERIAL  SCAN  ASSEMBLY  WITH  REFLECTIVE  DETECTOR  SURROUND  NARCISSUS 
AS  A FUNCTION  OF  OBJECT  POSITION  (SI)  AND  REAR  SURFACE  CURVATURE  (R) 

Negative  Case 


32 


Separation,  Si  - cm 


FIGURE  33 

SERIAL  SCAN  ASSEMBLY  WITH  REFLECTIVE  DETECTOR  SURROUND  NARCISSUS 
AS  A FUNCTION  OF  OBJECT  POSITION  (SI)  AND  REAR  SURFACE  CURVATURE  (R) 

Positive  Case 


acax  = 0.0087  cm  _____  r,  - *i  cm 

a cay  I .882  cm  r ^ * +°°cm 


FIGURE  34 

PARALLEL  SCAN  ASSEMBLY  WITH  REFLECTIVE  DETECTOR  SURROUND  NARCISSUS 

Posititive  Case 
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The  upper  curves  on  Figures  30  and  31  show  the  raster  ratio  which  was 
derived  for  an  infinitesimal  detector  dimension.  As  such,  it  gives  an  indi- 
cation of  the  spreading  or  expansion  of  the  detector  narcissus  image  with 
detector  field  lens  separation,  not  actual  narcissus  image  size.  The  data 
shown  was  expected  except  for  the  negative  curvature,  negative  separation  case 
shown  in  Figure  31.  Here,  an  inverted  image  occurs  at  the  detector  plane  at 
about  -1.8  cm  separation. 

Figures  32  and  33  show  the  narcissus  reduction  obtained  by  adding  a 
reflective  surround  to  the  detector  assembly.  This  reflective  material 
reflects  300°K  photons  from  the  detector  mechanical  assembly  into  the  relay 
optics,  thereby  making  the  detector  cold  assembly  appear  smaller.  Theoret- 
ically, this  technique  can  reduce  the  size  of  cold  area  to  that  of  the  detec- 
tor material.  As  a practical  value  in  this  case  however,  about  seven  times 
the  total  detector  array  area  is  assumed.  Results  using  this  scanner  show 
narcissus  is  insignificant  for  any  configuration. 

Figures  34  and  35  show  narcissus  for  the  parallel  scanner.  Here  narcis- 
sus is  acceptable  except  for  separations  of  zero  of  1 cm.  The  second  gener- 
ation detector  case  shown  on  Figures  36  and  37  shew  severe  narcissus  problems 
because  the  narcissus  returns  are  much  greater  than  the  NEP  unless  separa- 
tions are  maintained  over  2 cm. 

Based  on  the  above  data  a value  of  2 cm  wac  hosen  as  a desired  separation 
for  any  field  lenses  so  that  advanced  detectors  can  be  accommodated.  It  should 
be  emphasized,  however,  that  current  scanners  can  be  easily  accommodated  with- 
out any  significant  constraints  on  relay  design. 

NONLINEAR  LENS  NARCISSUS 

Narcissus  effects  due  to  elements  of  the  nonlinear  lens  were  calculated 
by  the  procedures  outlined  in  Appendix  C.  In  order  to  accomplish  this  task 
within  a reasonable  cost  and  time  schedule  an  early  lens  design  was  utilized 
which  is  slightly  different  than  the  lens  shown  in  Section  3.  In  addition, 
all  elements  are  considered  spherical  with  radii  equal  to  their  on  axis 
curvatures.  Neither  of  the  above  approximations  are  believed  to  introduce 
significant  error.  This  will  be  verified  during  early  analysis  during  the 
fabrication  phase. 

Results  of  the  narcissus  analysis  are  presented  on  Figure  38  for  each 
scanner  type.  They  were  obtained  by  scaling  the  detector  forward  to  the  non- 
linear lens  focal  plane  using  the  relay  parameters  developed  earlier. 

These  data  show  that  acceptable  narcissus  (significantly  belcw  NEP)  can 
be  obtained  for  the  current  technology  scanners  but  it  is  excessive  for  the 
second  generation  scanners.  This  may  be  due  to  the  fact  that  the  analytical 
technique  used  tends  to  be  a "worst  case"  for  most  elements  because  it 
assumes  that  all  of  the  energy  returning  to  the  narcissus  image,  that  lies 
within  the  on  axis  detector  ray  line,  will  illuminate  this  detector  (Figure 
C-5  of  Appendix  C) . For  elements  well  forward  of  the  image  plane  this  is 
quite  often  not  the  case.  Usually  much  of  the  energy  from  off-axis  detector 
cold  assembly  points  that  should  return  and  strike  the  detector  is  lost  on  the 
return  path.  This  is  illustrated  by  ray  trace  data  that  was  very  laboriously 
obtained  on  narcissus  from  the  extreme  front  surface.  This  narcissus  data 
is  listed  on  the  bottom  of  Figure  33  and  can  be  directly  compared  to  the  last 
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1 
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0.4536 
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Figure  in  parenthesis  is  ratio  of  narcissis  power  divided  by  NEP. 


FIGURE  38 

NARCISSUS  FROM  NONLINEAR  LENS  SURFACES 


value  of  the  approximate  data.  In  general,  narcissus  is  lower  by  more  than 
an  order  of  magnitude  when  computed  by  this  more  exact  technique.  An  excep- 
tion is  the  modified  serial  scanner  where  the  source  is  so  small  that  very 
little  energy  is  lost  on  the  return  path.  For  elements  located  close  to  the 
focal  plane,  the  probability  of  losing  energy  becomes  much  less  and  the 
approximate  method  for  computing  narcissus  (Appendix  C)  should  become  more 
accurate . 

NARCISSUS  CONCLUSIONS 

Narcissus,  while  being  somewhat  less  than  that  indicated  on  Figure  38, 
will  be  unacceptable  for  the  next  generation  FLIR  unless  dynamic  compensa- 
tion is  incorporated.  This  could  be  accomplished  on  a linear  scan  system 
by  storing  an  entire  scan  frame  and  subtracting  out  the  narcissus  bias. 
Another  technique  would  employ  a rotating  scan  array  to  take  advantage  of 
the  rotational  symmetry  inherent  in  the  narcissus  returns  of  the  nonlinear 
lens.  A full  two-dimensional  non-scanned  array  also  would  inherently  elimi- 
nate narcissus. 
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Section  3 

LENS  DESIGN  AND  OPTIMIZATION 

The  lens  design  phase  proceeded  along  lines  shown  in  Figure  39.  The 
effort  started  with  a materials  survey  to  find  the  best  material  for  con- 
struction of  the  infrared  nonlinear  lens.  Initially,  materials  suitable  for 
operation  in  either  the  3-5  micron  or  8-14  micron  regions  were  considered. 
Conclusions  were  to  use  silicon  for  the  3-5  micron  spectrum  and  germanium 
for  the  8-14  micron  spectrum  for  the  severe  aspherics.  The  remaining  elements 
could  be  made  of  other  materials  based  on  disperison  requirements.  The  pri- 
mary reason  for  the  aspheric  element  selection  was  advice  from  the  lens 
fabricator,  Frank  Cooke.  Both  materials  can  be  worked  as  easily  as  glass  and 
present  no  problems  of  fabrication.  In  addition  these  materials  can  be 
obtained  with  very  good  homogeneity  where  all  other  materials  are  questionable 
in  this  area.  Considering  the  thickness  of  the  large  aspheric  elements, 
homogeneity  is  of  upmost  importance.  Another  reason  for  selecting  silicon  and 
germanium,  are  their  low  dispersions  in  their  respective  operating  regions. 

The  dispersion  is  nearly  an  order  of  magnitude  smaller  than  that  of  any  other 
applicable  material.  Low  dispersion  greatly  simplifies  optical  design, 
especially  in  an  unconventional  design  such  as  this  one.  Finally,  both  mate- 
rials have  good  stable  mechanical  and  thermal  characteristics.  The  only  dis- 
advantage of  these  materials  is  their  high  index  of  refraction  but  this  does 
not  impose  any  restriction  on  design.  It  does,  however,  require  sophisticated 
optical  coatings  to  maintain  high  transmission.  The  required  coating  technol- 
ogy exists  and  does  not  outweigh  advantages  of  employing  the  high  index  mate- 
rials . 

As  the  analysis  progressed  the  8-14  micron  spectral  region  was  selected 
as  optimum  and  therefore  the  material  used  is  limited  to  germanium.  It  was 
later  shown  that  because  of  its  low  dispersion  the  entire  lens  could  be  made 
of  germanium. 

The  next  step  in  lens  design  was  to  initialize  the  existing  lens  computer 
program  with  parameters  of  the  IR  lens  and  establish  a first  order  design  for 
a single  wavelength.  An  index  of  4.003  was  used  with  an  8 inch  clear  aperture 
diameter  and  the  original  nonlinear  lens  distortion  function.  The  initial 
design  was  established  for  an  11  micron  wavelength  F/2.0  aperture. 

The  target  in  the  next  design  iteration,  was  a lens  with  an  optical  blur 
of  no  greater  than  16  microns  when  imaging  over  the  entire  8-14  micron  spec- 
tral region.  After  an  acceptable  design  was  achieved,  modifications  were 
attempted  to  produce  a well  behaved  optical  output  from  the  image  plane,  to 
facilitate  optical  relaying  to  the  scanner  and  detector.  In  general,  this 
consisted  of  establishing  a fixed  exit  pupil  for  the  F/2.0  ray  cone  output  by 
using  aspheric  field  lens  surfaces.  The  design  was  iterated  until  a fixed 
pupil  was  achieved. 

As  the  narcissus  analysis  progressed  in  the  Analytical  Phase  of  the  effort, 
it  became  apparent  that  narcissus  placed  severe  constraints  on  location  and 
curvatures  of  field  lens  surfaces.  Based  on  this  data,  it  was  necessary  to 
reiterate  the  design  while  forcing  the  field  lens  surfaces  away  from  the  lens 
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image  plane.  After  a successful  design  was  achieved,  it  was  noted  that  there 
was  a variation  in  exit  pupil  illumination  with  a change  in  field  angle.  This 
indicated  that  a non-uniform  F/number  existed  in  the  image  plane.  When  this 
problem  was  examined,  it  was  discovered  that  F/number  was  very  difficult  to 
define  on  the  nonlinear  lens  because  of  the  disparity  that  exists  between 
radial  and  tangential  focal  length.  A detailed  radiometric  analysis  was 
required  to  determine  how  to  define  the  lens  F/number.  This  analysis  is  pre- 
sented in  Appendix  D.  The  result  showed  that  the  best  way  to  define  F/number 
was  to  define  entrance  pupil  geometry  and  compute  an  equivalent  F/number  by: 

(24) 


where 

0 = Object  field  angle 

fr(6)  = Focal  length  in  the  radial  direction 

h(0)  = Radial  image  height 

a and  b = Major  and  minor  radii  of  the  entrance  pupil  ellipse. 

Using  this  equation,  the  magnitude  of  the  F/number  problem  was  established. 

It  was  found  to  be  excessive  and  led  to  reiteration  of  the  entire  design. 

Achieving  an  acceptable  F/number  uniformity  required  some  compromise  in 
relay  narcissus.  The  best  lens  design  achieved  during  this  contracted  effort 
is  shown  in  Figure  40.  Its  blur  and  F/number  performance  are  shown  on 
Figure  41.  Reiteration  and  further  optimization  of  thi6  design  are  continu- 
ing in  an  attempt  to  further  improve  F/number  uniformity.  We  believe  it  can 
be  made  acceptable  to  operate  with  current  generation  scanners,  if  not  by 
design,  then  by  designing  for  a smaller  F/number  and  stopping  down  within  the 
relay  to  achieve  uniformity.  As  a last  resort,  a radial  filter  or  variable 
thickness  element  coating  may  be  used  to  achieve  final  uniformity.  However, 
we  do  not  feel  this  will  be  necessary  with  current  generation  scanners  but 
will  definitely  be  required  for  second  generation  detectors. 

CONCLUSIONS 

The  infrared  nonlinear  lens  shown  in  Figure  40  is  feasible  and  can  fully 
support  the  VARVS  concept  with  either  current  or  second  generation  detector 
and  scanner  technology.  This  optimized  design  operates  in  the  8-14  micron 
spectral  region,  has  an  F/number  of  2.0  and  optical  quality  in  its  focal 
plane  is  represented  by  a one  sigma  blur  circle  radius  of  less  than  16  microns. 
System  performance  with  this  lens  is  represented  by  Curve  1 of  Figure  42  with 
a current  generation  serial  scanner  and  as  Curve  3 when  a second  generation 
monolithic  focal  plane  array  is  used  to  extract  the  image.  Actual  performance 
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may  be  somewhat  better  than  these  curves.  However,  because  the  blur  on  axis 
where  high  acuity  viewing  is  accomplished  is  less  than  the  16  micron  value, 
it  should  be  noted  that  the  16  micron  design  value  was  a "not  to  be  exceeded 
value"  anywhere  on  the  focal  plane.  Figure  41  showed  that  blur  on  axis  was 
actually  in  the  range  of  8-12  microns  with  an  average  of  10  microns.  The 
anticipated  performance  of  this  lens  would  be  as  shown  by  the  dashed  lines  on 
Figure  42.  The  result  is  significantly  better  thermal  resolution  at  all 
spatial  frequencies.  A comparison  of  thermal  resolution  at  the  human  visual 
dynamic  performance  spatial  resolution  point  (.86  £X£iS§.)  shows  the  performance 
improvements  listed  on  Figure  43.  The  improvement  is  very  substantial. 

The  analyses  indicate  narcissus  is  not  a problem  with  current  technology 
detectors  and  scanners.  However,  the  larger  cold  assembly  area  and  increased 
sensitivity  of  second  generator  detector  arrays  will  necessitate  some  method 
of  narcissus  compensation. 

Ho  fabrication  problems  are  anticipated  with  the  design  since  all  elements 
are  made  of  germanium  which  can  be  easily  ground  to  the  relatively  mild 
aspheric  shapes  required.  Use  of  germanium  also  minimizes  risk  in  both  fabri- 
cation, coating,  and  homogeneity.  Since  these  items  are  clearly  within  the 
state-of-the-art  which  is  not  true  for  many  candidate  materials. 

The  only  problem  with  the  design  is  the  maintenance  of  an  F/2.0  lens 
aperture  at  large  field  angles.  Improvement  made  in  the  last  few  weeks  of 
the  contracted  effort  convinces  us  that  further  refinement  in  design  can 
correct  this  problem.  Work  is  continuing  in  this  area  on  1RAD  funds. 
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MRT  PERFORMANCE  IMPROVEMENTS  AT  0.86  C/MR 
SPATIAL  FREQUENCY  BY  REDUCING  BLUR 
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Appendix  A 


BRIEF  DESCRIPTION  OF  THE  REMOTE  VIEWING  SYSTEM  (RVS) 

The  RVS  concept  is  based  on  the  fact  that  the  human  visual  capability  can  be 
represented  by  a resolution  capability  of  about  130,000  elements,  Drovided  that 
these  elements  are  sized  non-linearly  according  to  the  acuity  function  as  shown 
in  Figure  A-l.  An  image  with  this  characteristic  requires  only  about  2 MHz  video 
bandwidth  at  30  Hz  frame  rates.  In  comparison,  standard  techniques  would  require 
over  1,000  MHz  bandwidth  for  this  f ield-of-view  (180°)  and  resolution.  Even 
at  smaller  fields-of-view,  the  bandwidth  saving  is  significant.  A comparison  of 
bandwidth  requirements  for  varying  fields-of-view  for  the  conventional  linear 
acuity  function  and  for  the  RVS  foveal  concept  is  shown  in  Figure  A-2 . Appro- 
ximately two  orders  of  magnitude  decrease  in  BW  is  achieved  with  the  foveal  system 
at  FOV's  greater  than  20  degrees.  In  order  to  mechanize  the  concept  described 
above,  a method  must  be  devised  to  generate  an  image  which  satisfies  the  optical 
requirements  of  the  eye.  The  RVS  concept  contains  a lens  system  that  creates 
optical  "distortion"  by  varying  the  spacing  of  the  angular  resolution  elements 
to  duplicate  the  acuity  function  shown  in  Figure  A-l.  This  process  is  illustrated 
in  Figure  A-3.  The  lens  transfer  characteristic  required  and  the  technique 
for  reconstructing  the  image  at  a remote  location  is  also  shown  on  this  figure. 
System  operation  is  as  follows: 

The  image  transmission  system  scans  the  photocathode  of  the  vidicon  or  photo- 
detectors of  an  imaging  array,  transmits  this  signal  to  the  remote  location,  and 
recreates  the  image  on  a CRT  or  light  valve  tube.  In  the  original  RVS  concept, 
the  distorted  image  is  expanded  using  a lens  system  with  a transfer  characteri- 
stic identical  to  the  sensor  lens  and  imaged  on  a spherical  screen  concentric 
with  the  nodal  point  of  the  lens. 

Obviously,  for  the  above  image  transmission  system  to  perform  adequately,  the 
optical  axes  of  both  the  sensor  and  projector  must  have  the  same  alignment  as 
the  viewer's  eye.  The  initial  RVS  system  concept  used  the  approach  outlined  in 
Figure  A-4.  The  position  of  the  projector  is  slaved  to  the  camera  by  a high 
accuracy  position  servo,  with  the  camera's  angular  position  commanding  the  pro- 
jector's position  relative  to  fixed  ground  station  reference  coordinates.  The 
viewer  at  the  ground  station  thus  has  the  same  angular  perspective  as  he  would 
if  he  were  located  in  the  remote  vehicle.  The  sensor  and  projector  must  also  be 
aligned  with  the  viewer's  foveal  axis.  In  the  original  concept  a Honeywell 
oculometer  was  employed  for  this  function.  The  oculometer  measures  the  angle 
between  the  eye's  foveal  axis  and  the  projector's  optical  axis.  This  error  signal 
is  transmitted  to  the  remote  vehicle  and  commands  the  camera  to  move  until  the 
angular  error  is  reduced  to  zero.  As  the  camera  moves,  the  projector  follows 
through  the  slaving  loop.  The  control  mode,  presently  under  study,  is  somewhat 
different,  however.  The  observer's  head  position  instead  of  his  eye  position  is 
utilized  to  point  the  remote  camera.  The  operational  difference  resulting  from 
this  simplification  is  that  when  the  viewer  uses  his  peripheral  vision,  he  must 
learn  to  rotate  his  head  towards  the  area  of  interest  rather  than  his  eyes.  A 
reticle  may  be  required  to  show  the  observer  the  location  of  the  highest  acuity 
area  of  the  display. 
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Figure  A-1.  Human  Eye  Characteristics 
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Appendix  B 


APPLICATION  OF  THE  NIGHT  VISION  LABORATORY  (NVL) 
THERMAL  VIEWING  SYSTEM  STATIC  PERFORMANCE  MODEL  TO 

THE  RVS 


It  was  suggested  that  the  NVL  Thermal  Viewing  System  Static  Performance  Model 
be  used  to  evaluate  the  performance  of  the  Remote  Viewing  System  (RVS).  However, 
repeated  attempts  to  convert  the  RVS  parameters  directly  to  the  NVL  model  have  led 
to  the  following  problem.  The  radial  distortion  function  of  the  foveal  lens  does 
not  lend  itself  to  an  MTF  analysis  as  a function  of  object  field  angular  spatial 
frequency  as  called  for  in  the  NVL  model.  All  parameters  can  be  converted 
successfully  except  for  the  scan  velocity  term  because  a linear  raster  scan  on 
the  lens  image  plane  will  create  a variable  angular  velocity  and  variable  direction 
scan  in  the  object  field.  This  is  depicted  in  Figure  B-l.  Extreme  complexity 
results  when  attempts  are  made  to  convert  spatial  into  temporal  frequency.  This 
is  illustrated  by  the  rotation  of  the  f bar  pattern  in  the  lens  image  plane 
shown  in  Figure  B-l.  Given  enough  time,  an  analysis  could  be  made  in  a manner 
compatible  with  the  NVL  model.  However,  the  analysis  is  much  simpler  if  performed, 
not  in  object  field  angular  frequency  (cycles /milliradian)  but  in  spatial  frequency 
terms  (lines/millimeter).  For  our  purpose  of  optimizing  the  RVS  lens,  it  is 
simpler  to  work  in  terms  of  spatial  frequency  on  the  foveal  lens  focal  plane. 
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This  simplicity  arises  because  seven  of  the  nine  MTF's  are  independent  of 
object  field  angle  at  this  foveal  lens  focal  plane  location,  and  the  scan 
velocity  is  undirectional  and  uniform  at  this  location,  thereby  making  easy 
conversions  from  spatial  to  temporal  parameters.  The  only  non-linear  conversions 
necessary  are  simple  geometrical  ones  which  translate  from  focal  plane  to  object 
field  and  display  space.  The  advantages  of  working  in  the  spatial  frequency 
terms  will  become  clear  as  the  analysis  is  developed.  In  the  following  develop- 
ment, the  NVL  model  approach  will  be  used  precisely  but  will  be  applied  in  the 
foveal  lens  focal  plane  as  a function  of  linear  spatial  frequencv  (1/mm).  Parameters 
will  be  covered  in  the  same  order  as  they  are  in  the  NVL  Report  ^ , which 

describes  the  model  in  detail. 


B.l  MTF's 

1.  Optical  MTF  The  optical  MTF's  consist  of  a diffraction  MTF  and  a Gaussian 
MTF. 

(a)  Diffraction  In  angular  terms,  the  diffraction  MTF  is  referenced  as  Equations 
(9)  and  (10)  of  the  NVL  report: 


H _(f  ,9)  “ — fcos  1 A - A(1  - A^)^^l 

opt  x tt  *-  J 

where  A - 


(B-l) 

(B-2) 


where  L(6)  is  the  equivalent  focal  length  which  changes  over  a 50/1  range  as 
object  field  angle  6 changes.  The  angle  6 is  the  absolute  angle  between  the 
point  of  interest  and  the  lens  optical  axis.  At  the  foveal  lens  image  plane 


f 


(B-3) 
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(b)  Blur  - A similar  simplicity  exists  here.  The  MTF  equation  with  the  angular 
term  b of  Equation  (11)  of  Reference (B-l) replaced  with  its  equivalent  is: 


**blur 


6) 


, 2 2 
uTT  0 


L(9)‘ 


(B-7) 


The  foveal  lens  inherently  has  a constant  spatial  blur  over  its  entire  focal  plane, 
so  that  the  sigma  (a)  of  Equation(B-7) is  a constant.  Substituting  Equation (B-5) into 
(B-7)  we  see  the  blur  MTF  simplifies  to 

Hblur(Sx)  = 6XP  C-  21,2  (B-8> 

Thus  this  MTF  like  the  diffraction  MTF,  is  no  longer  a function  of  object  field 
angle  because  the  focal  length  variable  has  been  removed. 

2.  Detection  MTF  - The  spatial  filter  MTF  of  the  detector  is  defined  as: 


Sin(irf  Ax) 

(f  ,0)  = ; . X — = Sinc(f  Ax) 

Det  x’  irf  Ax  x (B- 

x 

It  is  also  complex  in  our  system  because  the  angular  projection  of  the  detector 
into  the  object  field  (Ax)  in  this  equation  varies  with  absolute  object  field 
angle  (8).  Since  the  detector  height  is  still  uniform  at  the  lens  focal  plane, 
shown  in  Figure (B-2) as  Ah,  Equation(B-9) can  be  restated  as: 


Sin(irS  Ah  ) 

R (S  ) = *■  -X-  (B-10) 

Det  x ttS  Ah 

X x 

Again  the  MTF  becomes  independent  of  object  field  angle.  Note  from  Figure(B-2) that 
the  detector  height  (Ah^)  is  a function  of  detector  size(a),  detector  system  focal 

length  (L  ) , and  relay  focal  length  (L  ) , viz: 

D ^ 


Ah  = a ^ (B-ll) 

x x 


There  are  two  ways  to  arrive  at  the  required  and  LQ  based  upon  whether  the  scanner 

or  detector  characteristics  are  known.  If  the  scanner  is  available,  horizontal  FOV 

(0  ) and  instantaneous  FOV  (a)  are  known.  Then  the  detector  focal  length  is 
s 
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FIGURE  B-1  GP7.  1037  .09 

SCAN  DISTORTION  INTRODUCED  BY  FOVEAL  LENS 

where  is  the  image  plane  spatial  frequency  and  f is  its  object  field 
angular  equivalent  measured  along  the  scan  line  projection  in  the  object 
field  (y  direction  on  Figure  B-1).  Solving  for  f in  Equation (B-3)and 
substituting  this  for  f^  in  Equation  (B-2).  y 

A = X F„  S (B-4) 

Since  the  F/number  of  our  lens  is  constant,  the  diffraction  MTF  is  no  longer  a 
function  of  object  field  angle.  Thus  we  may  write  H (S  ) which  indicates 
that  the  MTF  is  a function  of  the  independent  variabfl  S Xonly.  Note,  however, 
that  conversion  to  object  field  angular  spatial  frequency  is  very  simple 
because  focal  length  is  constant  over  small  angular  increments  and  may  be 
determined  from 

fy  - SxL(6)  (B-5) 

where  y is  along  the  scan  line  projection  in  the  object  field 
likewise 

f - SL(8)  (B-6) 

w y 

where  w is  normal  to  the  scan  direction  in  the  object  field 
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FIGURE  B-2 

SCANNER  RESOLUTION  ELEMENT  GEOMETRY 

On  the  forward  end  of  the  scanner  we  must  cover  the  foveal  lens  image  plane 
height  of  2h(0  ^)  which  subtends  the  angle  6g.  Thus, 


Solving  for  Lc 


h(6  ) _ 9 

max  tan  s 


h(0  ) 

max 

tan  6 /2 
s 


(B-13) 


(B-14) 


If  the  detector  characteristics  are  known,  the  focal  lengths  are  a function 
of  detector  size  (Ah)  projected  unto  the  image  plane  as  shown  in  Figure  (B-2). 
Detector  size  Ah  can  be  computed  directly  from  either  the  on-axis  resolution 
required,  the  number  of  scan  lines  required  across  the  vertical  FOV,  or  bandwidth/ 
response  restrictions  and  frame  rate  requirements.  The  focal  lengths,  and  L^, 

are  then  selected  to  make  the  detector  dimension  appear  as  the  required  Ah  at 
the  foveal  lens  focal  plane.  The  scanner  horizontal  FOV  given  by  scanner  mechanics 
must  cover  the  total  image  plane  height  (2h(0max)).  1°  either  case  the  detector 
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MTF  becomes : 


>W(SX>  - Sine  (B-15) 

x D 

Again  this  >OF  is  independent  of  object  field  angle. 

3.  Detector  Electronics  MTF  - It  is  in  the  MTF,  the  detector  electrical 
response,  that  we  get  into  real  trouble  trying  to  work,  in  object  field 
angular  space.  For  a conventional  linear  optical  system,  a linear  detector 
scan  velocity  converts  into  a scaled  but  linear  angular  scan  in  the  object 
field.  This  is  not  true  in  our  system  as  was  shown  in  Figure  B-l.  A linear 
scan  in  the  x direction  on  the  image  plane  results  in  angular  velocities  in 
both  0x  and  0^  directions  in  the  angular  object  field.  Both  of  these  angular 

components  are  nonlinear  functions  of  both  x and  y position  on  the  image 
plane.  Thus,  converting  from  spatial  frequency  to  temporal  frequency  becomes 
very  complex.  All  of  this  can  be  avoided  by  working  in  linear  spatial  plane 
terms.  If  the  scanner  has  an  angular  scan  velocity  6,  then  the  linear  motion 
of  the  instantaneous  FOV  on  the  foveal  lens  image  is 

V = 8L  (B-16) 

x C 

The  conversion  to  temporal  frequency  (f)  is  therefore 

f = V S (B-17) 

x x 

This  is  a constant  conversion  and  not  a function  of  time.  Therefore,  all 
electronic  MTF’s  of  the  NVL  model  are  valid.  These  are 
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4.  Display  - The  RVS  display  is  the  inverse  of  the  foveal  lens,  which  results 
in  a conventional  linear  raster  generated  on  the  CRT.  The  CRT  has  a constant 
spot  size  and  the  expansion  optics  has  a constant  blur  at  the  object  focal  plane. 
Again  this  MTF,  if  derived  in  the  linear  spatial  plane,  will  not  be  a function 
of  object  angle.  If  the  optical  blur  and  CRT  spot  size  are  combined  and 
assumed  to  have  a Gaussian  MTF,  a composite  sigma  (o^)  results  and  the  MTF  is: 


HDisp(Sx)  = eXp  l>2ir2(rod)2sx2  ] 
where  r is  the  physical  ratio  of  format  sizes  ; viz 


r 


HlENS  image 
Display  crt 


(B-18) 


(B-19) 
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By  contrast,  if  this  were  accomplished  in  the  object  angular  plane,  the  MTF 
would  be  much  more  complex,  viz 


2ir2(rod) 

L(0)2M 

where  M is  any  system  angular  magnification  from  object  field  to  the  viewer. 
Again  the  simplicity  is  obvious. 


**Disp  (fx,6  ,M^ 
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(B-20) 


5 & 6.  Stabilization  and  Eyeball  - The  remaining  two  MTF's  are  the  only 
two  that  are  not  simplified  by  working  in  linear  spatial  rather  than  angular 
terms.  First,  stabilization  tends  to  be  angular  input  to  the  system.  Using 
the  MTF  from  the  NVL  report: 


4 


■Wf*>  • “■>  <-«„  > 

Converting  to  the  foveal  lens  image  plane  results  in 

■WV’  - «p  [-ps,2^*)2] 

Similarly,  the  eye  views  the  display  in  angular  terms.  The  NVL  MTF  is 


HEye(fx)  = expj 

Equation(B-23)must  be  converted  to  the  foveal  lens  image  plane 


HEye(Sx,6)  = 6XP 


al  lei 

[-3T] 


(B-21) 


(B-22) 


(B-23) 


(B-24) 


In  conclusion,  seven  MTF’s  have  been  simplified  at  the  expense  of  two  that 
have  been  made  slightly  more  complex  by  the  conversion  to  linear  spatial 
frequency. 


B.2  NEAT 

NEAT  can  be  used  as  defined  in  the  NVL  document  since  it  is  independent  of 

spatial  frequency,  therefore  no  changes  are  required.  However,  it  is  a function  of  F/ 

number  and  a discussion  is  in  order  concerning  which  is  the  correct  F /number  i 

to  use  in  the  calculation.  It  is  not  simply  the  foveal  lens  F /number  or 

objective  lens  F /number  as  stated  in  the  NVL  document.  For  NEAT  calculations,  i 

the  F /number  must  represent  the  actual  ray  cone  supplying  the  detector. 

In  general  our  lens  will  be  limited  to  a fixed  image  height  h(9  ) 

max 

that  is  based  on  fabrication  considerations.  The  scanner  must  cover  this 
total  image  in  the  most  efficient  manner.  The  optical  geometry  of  the  lens 
and  scanner  can  be  reduced  to  the  basic  arrangement  shown  in  Figure (B-3). 

On  this  figure  the  ray  cone  that  supplies  the  detector  is  Note  that 

a maximum  ray  cone  (d>^)  that  the  detector  can  utilize  exists.  Obviously, 


■v  „ 


i 
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the  best  design  would  have  This  is  not  always  possible,  however, 

because  the  lens  may  be  required  to  operate  with  available  scanners.  Therefore, 
it  will  be  necessary  to  define  both  of  these  ray  cones  and  utilize  the  smallest 
for  NEAT  and  diffraction  MTF  calculations.  Before  processing,  note  that  ray 
cones  may  also  be  defined  by  F/numbers  as  shown  in  Figure (B-4). 


D <>> 


Figure  B-4.  Ray  Cone  Parameters 
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Since  anywhere  along  the  ray  bundle 


F#  = L/D 

<t>  = 2 tan-1(D/2L) 


(B-25) 


and  F/number  may  be  defined  by 


F(/  “ 2 tan  <j>/2  (B-26) 

Therefore,  we  can  use  F/number  instead  of  $ to  define  ray  cones.  By  the  nature 
of  collimation,  the  ray  bundles  between  the  collimating  lens  and  the  detector 
will  be  of  constant  cross  section,  i.e., 


And 


Therefore 


And 


F#1  F#  3 

F - f 
*#4  “ Lc  F#2 


(B-27) 

(B-28) 


The  larger  of  these  two  F/ numbers  must  be  used  for  both  diffraction  and  NEAT 
calculations.  Now  note  that  F^  is  the  F/number  of  the  foveal  lens  while  F^ 

is  tha_t  generated  by  the  scanner  internal  aperture  stop  (Dg)  and  we  may  write: 


(B-29) 


For  any  particular  scanner,  Dg  will  be  known.  From  Equation (B-28) the  smallest 

F/number  which  translates  into  the  largest  ray  cone  we  can  supply  the  detector 
with  is: 


which  is  also  obvious  from  inspection. 


(B-30) 


i 
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A 


For  our  ONR  work,  a scanner  Has  been  specified  while  for  the  USAF  effort  only 
a detector  geometry  ia  available.  The  following  derivations  will  try  to  accommodate 
both  arrangements.  If  the  scanner  is  specified,  Dg  and  will  be  known.  If 

is  not  given  it  can  be  calculated  from  the  scanner  instantaneous  FOV  (a)  and 
detector  size  (a)  in  the  same  plane.  Then 


L 


D 


(B-31) 


Also,  if  the  scanner  is  specified,  its  maximum  scan  angle  (0  ) shown  in  Figure 

(B-2).is  known.  Since  the  scanner  must  scan  the  entire  image,  6^  must  cover 

2h(0  ).  Then 

max 


L 


C 


h(6  ) 

max 


(B-32) 


It  is  now  possible  to  compute  F,,,  from  Equation  (B-27) 


#4 


#4 

a tan  9 /2 

-Z § — 

a h(0  ) 

y max 


= F 


#2 


(B-33) 


where  F„  is  the  foveal  lens  F/number.  Remember,  for  the  diffraction  MTF  and 
NEAT  calculations,  we  use  the  larger  of  F^  or  F^. 

If  a scanner  is  not  defined,  one.  must  be  theorized.  This  can  easily  be 
accomplished  by  procedures  outlined  in  Lloyd  B-l.  Generally,  scanning  hardware 
(mirror,  drum,  etc.)  are  selected  to  accomplish  a given  scan  angle  6g  (number  of 

facets).  Collimation  lens  focal  length  (L^,)  can  then  be  computed  by  Equation  (B-32). 

From  the  required  scan  lines  across  the  maximum  image  height  h(®max)>  the  apparent 

detector  size  (Ah)  at  the  image  can  be  calculated.  The  angular  instantaneous  FOV 
of  the  scanner  can  be  calculated  as  follows,  using  the  ori-axis  geometry  of  Figure 

(b-2): 


a 

y 


Lp  can  then  be  calculated  by 


L 


D 


a* 


y 


D 

s 


can  then  be  defined  by  matching  the  lens  F/number 


(B-34) 


(B-35) 


(B-36) 
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B.  3 MKT  CALCULATIONS 


The  folic  ing  MRT  equation  modifications  are  required  so  that  the  computation 

may  be  performed  in  linear  spatial  frequency  terms.  First,  in  the  NVL  MRT  equation,  Ay 

must  be  replaced  by  the  apparent  detector  size  at  the  foveal  lens  image  plane, 

i.e.,  it  must  be  the  Ah  defined  on  Figure (B— 2).  As  previously  demonstrated  in 

Equation  (B-ll), 


(B-37) 


Also,  in  the  MRT  equation,  it  is  best  to  compute  the  Q integral  in  terms  of 
temporal  frequency.  This  eliminates  the  velocity  term  in  the  MRT  equation 
and  makes  the  Q integral  easier  to  compute.  The  Q integral  is  therefore 


Q<£'6)  ■/’  Hfy  <B~38) 

Of  these  terms,  only  H^,  the  transfer  function  for  a rectangular  bar  of  width  w, 

has  not  been  defined.  This  transfer  function  is  in  linear  rather  than  angular 
dimensions,  i.e., 

Sinc^w£xj  = Sine  ^WSxj  (B-39) 

where  ^ 

w = 2i—  (B-40) 

x 

The  MRT  equation  written  to  show  the  dependency  of  two  variables  is 


MRT(Sx,0) 


SNRtt2NEAT 

MTFrom(!Ve> 


Ah  S Q(f , 0) 

y x * 

AfNFR  Ce  nOVSC  J 


1/2 

(B-41) 


This  equation  results  in  an  MRT  very  weakly  dependent  on  8.  To  obtain  the  MRT 
for  any  field  angle  6 we  convert  the  spatial  frequency  term  into  an  angular 

frequency  term  by  using  Equation (B-5)containing  the  focal  length  function: 


f 

W 


= SxL(0) 


(B-5) 


Note  this  will  be  the  angular  spatial  frequency  in  the  scan 
bars  normal  to  the  scan  direction).  It  could  be  related  to 

this  does  not  appear  to  be  required  at  this  point. 


direction  (target 

f and  f but 
x y 
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Appendix  C 


THEORY  OF  NARCISSUS  ANALYSES 

The  problem  of  the  IR  detector  "seeing"  cold  focal  plane  surfaces  which  are 
reflected  by  the  system  optics  in  addition  to  the  warm  target  is  labeled  the 
"Narcissus  Problem".  The  source  of  the  narcissus  problem  is  illustrated  in  Figure 
C-l.  In  Figure  C-l  the  detector  assembly  is  shown  at  the  center  scan  position. 

All  areas  of  the  rear  clear  aperture  of  the  optics  are  flooded  with  photons  emitted 
from  the  300°K  structure  surrounding  the  detector  assembly.  Most  of  these  photons 
are  absorbed  by  the  lens  structure  or  pass  out  into  the  object  field.  A small 
number  however,  are  reflected  back  to  the  image  plane  by  the  various  lens  elements 
and  therefore  some  of  these  returning  photons  fall  on  the  detector  to  alter  its 
electrical  output.  Alternatively,  the  emitted  photon  flux  field  from  the  warm 
structure  may  be  thought  of  as  being  disrupted  by  the  cold  detector  assembly  which 
for  all  practical  purposes  has  no  emission.  Therefore  a cold  "hole"  is  introduced 
into  the  flux  field.  The  narcissus  problem  is  generated  when  this  hole  is  made 
to  translate  across  the  optical  image  as  illustrated  by  the  dashed  lines  in  Figure 
C-l.  This  occurs  during  the  mechanical  scanning.  This  translating  region  of  non- 
emission will  cause  a fluctuation  in  the  reflected  photon  distribution  at  the 
detector  plane.  If  this  spatial  fluctuation  is  smaller  than  the  total  scan,  an 
unwanted  signal  will  be  generated  during  the  scan.  The  generation  and  prediction 
of  the  narcissus  image  are  described  below. 

The  best  way  to  compute  significance  of  the  "narcissus"  energy  is  to  assume 
that  the  detector  cold  assembly  is  emitting  at  300°K  and  trace  the  reflections  that 
return  to  the  detector  plane.  The  spatial  distribution  of  reflected  portion 
of  this  energy  in  the  detector  plane  is  the  photon  energy  that  will  be  absent 
in  the  actual  situation,  hence,  a "cold"  output  will  be  obtained  from  the  detector. 
This  process  is  shown  in  Figure  C-2. 

At  the  "start  scan"  position  shown  in  Figure  C-2 (a),  three  points  of  emission 
on  the  detector  cold  assembly  are  shown.  In  the  most  severe  situation  (negative 
magnification)  the  reflected  energy  will  appear  on  the  opposite  end  of  the  scan. 

For  example,  if  the  optics  are  theorized  to  form  an  image  somewhat  in  front  of 
the  focal  plane  by  reflection,  the  three  ray  bundles  will  converge  forward  of  the 
focal  plane  and  expand  towards  the  focal  plane.  This  energy  forms  a "blur"  area 
in  the  focal  plane  as  shown  in  Figure  C-2 (a).  Remember  this  area  represents  a 
region  of  "less  photons"  than  the  surrounding  area.  Since  the  detector  does  not 
lie  within  this  area,  its  output  is  not  influenced  by  it. 

In  Figure  C-2(b)  the  detector  cold  assembly  is  shown  at  1/4  scan.  Note 
that  here  the  "narcissus"  image  has  moved  towards  the  detector  but  still  has  not 
influenced  its  output.  In  Figure  C-2(c)  the  detector  is  shown  at  center  scan. 

Here  it  lies  directly  within  the  "narcissus"  area  and  will,  therefore,  have  less 
output  or  will  have  a cold  output  created  by  the  lower  photon  input.  As  the 
detector  completes  the  scan,  it  moves  outside  of  the  "narcissus"  area  as  shown 
in  Figure  C-2(d)  and  is  no  longer  affected  by  narcissus.  The  total  result  then 
is  a "cold"  image  resembling  the  narcissus  blur  area  at  the  center  of  the  FLIR 
display.  This  is  illustrated  in  Figure  C-3. 

To  compute  the  true  effect  of  narcissus,  a ray  trace  of  all  points  within 
the  detector  cold  assembly  should  be  conducted.  Sufficient  ray  tracing  must  be 
conducted  to  map  the  power  density  in  the  detector  plane.  Much  of  this  laborious 
procedure  can  be  avoided  by  a series  of  simplications . For  example,  if  the  portion 
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Figure  C-1.  Source  of  Narcissus  Energy 


'For  analysis  purposes  only.  The  Cold  Detector  does  not  emit.  The  background  would  emit  this 
energy  if  detector  assembly  was  absent. 

GP77  0347  4 


FIGURE  C-2 

NARCISSUS  IMAGE  GENERATION 


6< 


FIGURE  C-3 

NARCISSUS  DISPLAY  OF  FIGURE  C-2 

of  the  narcissus  image  that  contributes  to  illumination  of  the  on-axis  detector  can 
be  established,  simple  optical  and  radiometric  theory  can  be  utilized  to  define  the 
maximum  strength  of  the  narcissus  effect.  A similar  procedure  can  be  used  to 
estimate  the  size  of  the  narcissus  image.  This  is  developed  in  the  following 
paragraphs. 

On  Figure  C-4(a),  emission  of  the  on-axis  point  of  the  cold  assembly  is 
illustrated.  If  narcissus  analysis  is  being  performed  on  a surface  forward  of 
the  scan  mechanism,  the  only  energy  that  will  be  reflected  back  to  the  detector 
is  confined  within  the  F/number  cone  of  detector.  This  cone  is  defined  by  the 
basic  optical  design.  From  elementary  geometric  optics,  the  energy  returning 
from  the  forward  surface  will  be  confined  to  this  cone  but  altered  by  the  optical 
magnification.  With  reference  to  Figure  C-4  it  is  easily  seen  that  the  blur 
circle  radius  is: 


R = S 


S 

2 M FNo 


(C-l) 


where  S = Distance  narcissus  image  is  in  front  of  detector  plane 
M = Magnification 

The  above  are  shown  in  Figure  C-4 (a)  and  (b)  and  is  the  narcissus  power  contribution 
for  a single  point  on  the  cold  emitting  assembly.  Obviously,  the  detector  will  be 
affected  by  the  integrated  results  of  the  emission  of  all  points  on  the  cold  assembly 
that  can  pass  through  the  optical  system  and  be  reflected  back  to  the  detector. 

Three  of  these  are  shown  in  Figure  C-4(c).  It  is  easily  seen  that  all  of  these 
points  contribute  to  output  of  the  detector  located  near  the  optical  axis.  The  difficult 
task  is  to  estimate  just  how  much  of  narcissus  image  is  capable  of  illuminating  the 
detectors.  This  is  a function  of  direction  of  the  chief  rays  forming  the  narcissus 
image . 


In  Figure  C-5  the  image  cone  of  the  on- axis  detector  is  shown.  For  a correctly 
designed  optical  system,  it  will  be  impossible  for  any  returning  ray  to  reach  the 
on-axis  detector  unless  it  falls  within  this  cone.  Therefore,  rays  forming  the 
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NARCISSUS  BLUR  CIRCLE  RADIUS 


N = 300° K Radiance  of  Detector  Cold  Assembly 
t - Optical  Transmission  Reflectance  Product 
rN  = Radiance  of  Narcissus  Image 
Wq  = Detector  Area/S2  = AD/S2 
PD  = rNA;  Wq 

PD  = 4 t NADXY 

S2 

Narcissus  Image  Area  that 
can  Illuminate  an  On-Axis 
acay  Detector  A|  = 4 x y 

y.  In  Most  Cases 

rNAD  M acax 


PD=- 


S FNo 


or  M Whichever  is  Smaller 

2 FNo  2 


OP77  0M7  7 

FIGURE  C-5 

ON  AXIS  NARCISSUS  POWER  CALCULATION 


62 


narcissus  image  cannot  illuminate  this  detector  if  they  occur  farther  from  the  optical 
axis  than  the  distance  y on  Figure  C-5.  The  narcissus  area (A.)  on  Figure  C-5 
is  therefore  1 


= 4 x y 


(C-2 ) 


Now  employing  another  simple  optical  theory,  that  the  image  radiance  (Nj_)  is  numer- 
ically equal  to  the  object  brightness  (N)  times  an  attenuation  factor  T.  which  is 
the  product  of  all  transmissions  and  reflectances  throughout  the  narcissus  path.  The 
derivation  of  the  on-axis  power  equaiton  is  shown  on  Figure  C-5.  For  a linear 
array,  the  last  equation  almost  always  applies,  i.e.: 

T N A_  M a 
r,  D cax 

PD 2 (C-3) 


Unfortunately  on-axis  power  is  not  sufficient  to  fully  assess  the  significance 
of  narcissus.  The  spatial  structure  on  the  display  must  be  known,  because  in 
many  cases  the  narcissus  power  may  be  many  times  the  detector  NEP  but  its  influence 
may  cover  all  detectors  equally  during  scan.  Therefore  it  will  not  be  seen.  This 
is  especially  possible  if  positive  magnification  exists. 


Starting  with  the  more  simple  case  of  negative  magnification,  the  case  depicted 
on  Figure  C-2,  estimation  of  narcissus  image  size  is  very  easy.  For  S > a , 
the  blur  image  width  (w)  in  the  detector  plane  is  cax 


w 


2R 


S 

M FNo 


(C-4) 


The  velocity  of  this  blur  during  scan  is 


blur 


- MV. 


(C-5) 


Since  the  detector  moves  in  the  opposite  direction  with  velocity  V,  the  relative 
velocity  between  the  two  is 

Vf  = V (1  - M)  (C-6) 

The  time  required  to  scan  the  narcissus  image  is 
= 2R  = 2R 

Slur  = V(l-M)  (c-7> 


The  time  required  to  scan  the  total  format  is 


t 

s 


X 

V 


(C-8) 


where  X is  format  width 


The  ratio  of  the  display  width  covered  by  the  narcissus  image  is 

Slur  _ 2R 

r t (l-M)x 

s 


(C-9) 


For  positive  magnification,  the  theory  is  much  more  complex.  First  it  is 
necessary  to  know  the  nature  of  the  basic  image  to  proceed  with  this  calculation. 
For  a telecentric  image  (all  chief  rays  parallel  to  the  optical  axis).  The  theory 
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can  be  developed  as  follows.  Figure  C-6  must  be  used  for  this  derivation.  This 
figure  is  a view  in  the  x plane.  Therefore  the  detector  and  Narcissus  image  are 
vertical  lines  normal  to  the  paper  that  move  from  the  bottom  of  the  page  (start 
scan)  to  the  top  (end  scan).  In  this  case  (positive  magnification),  the  narcissus 
image  goes  in  the  same  direction  as  the  detector,  and  may  continue  to  illuminate 
all  detectors  throughout  the  scan,  unless  vignetting  occurs  within  the  optical 
system.  Again,  by  assuming  a telecentrlc  optical  image,  an  assessment  can  be 
made  of  this  effect  on  the  displayed  image.  In  this  case  the  optics  will  limit 
narcissus  rays  to  the  format  diagonal  dimension  on  the  detector  plane  because  the 
optics  were  designed  to  do  this.  This  means  that  no  ray  can  rise  higher  than  the 
height  (W)  in  the  narcissus  image  plane,  again  remembering  that  moving  vertical 
pencils  are  being  considered.  In  this  case,  the  narcissus  energy  will  be  completely 
vignetted  when  the  detector  array  is  at  a height  (— ) from  the  optical  axis.  It 
seems  logical  therefore,  to  assume  the  image  disappears  at  this  point.  The 
derivation  shows  that  for  this  assumption  the  ratio  of  narcissus  occupied  display 
format  height  is: 


r (S) 


1.41  + S 

M M a FNo 
cay 


(C-10) 


A similar  analysis  for  negative  values  of  S yields 


r (-S) 


1.41 

M 


S 


M a FNo 
cay 


Image 

Diagonal 
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FIGURE  C-6 

NARCISSUS  SIZE  FOR  POSITIVE  MAGNIFICATION 
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Appendix  D 

NONLINEAR  LENS  F/NUMBER  DETERMINATION 


The  derivation  of  the  F/number  of  the  nonlinear  lens  is  defined  in  this 
Appendix.  Germaine  terms  are  defined  on  Figure  D-l.  Photons  are  radiated  from 
the  earth  (or  target)  at  a slant  distance  (S)  and  radial  distance  (P.)  from  the 
nonlinear  lens.  The  variables  S and  R are  much  greater  than  the  clear  aperture 
diameter  (DCA) . This  photon  flux  is  specified  by  a radiance  (N)  in  watts/CM  -Steradian. 
For  an  infinitesimal  area  in  the  object  field  (A^),  the  power  (Pa)  entering  the 
lens  pupil  of  area  (A)  is: 


where 


P = NA  W 
a t 

id  = Solid  angle  subtended  by  the  entrance  pupil 


NA  A 
t 


This  entire  flux  must  appear  on  the  lens  image  plane  within  the  image  (A^)  of 
object  area  (A^) . 

A^  = dt  dr 

The  power  density  of  the  focal  plane  is  therefore 
P NA  A 

V = -A  - ~r 

i S dt  dr 

Now  A can  be  related  to  the  image  plane  through  lens  focal  length.  By 
definitionthe  radial  focal  length  (f^)  is 

1 r<»)  ' 3?  <“-5> 

and  the  tangential  focal  length  (f^)  is 

ft<6)  - f£-  (D-6) 

t 

where  0^  and  0^  are  the  radial  and  tangential  angular  subtense  of  At.  Therefore, 

dr  - fr(9)d(er)  (D-7) 


dt  - ft(9)d(0t) 


Now  A, 


S d0  d0 
r t 


Substituting  this  into  Equation  (D-4)  we  have 


W = 


N A d6  d6 
r t 

dt  dr 


Now  substituting  Equation  (D-8)  and  (D-9)  into  (D-10) 


W = 


N A 

fr(6)ft(0) 


(D-10) 


(D-ll) 


The  focal  length  fr(0)  is  available  from  the  original  image  height  function  used 
in  lens  design  in  the  form  of  h(9)  vs  0,  thus 


£r<e>  - ^ 


(D-12 ) 


The  tangential  focal  length  is  not  so  obvious.  It  must  be  calculated  by 
circular  symmetry  considerations.  The  same  number  of  area  elements  must. exist 
around  the  object  field  circumference  2irR  as  exist  around  the  image  plane  cir- 
cumference 2tt  h(0).  In  the  object  field  the  number  of  elements  is: 


N 

eo 


2ttR  _ 2 tt  Sin  9 

d0  S d9 

t t 


(D-13) 


In  the  image  plane  the  number  of  elements  is 

N . -***§! 

ei  dt 


Then  Equation  (D-13)  and  (D-14)  may  be  set  equal  to  each  other 


Sin9  _ h(9) 
d0t  " dt 


Or  by  definition  of  the  focal  length 


ft(9) 


dt  _ h(9) 
d0t  ~ Sin0 


(D-14) 


(D-15 ) 


(D-16) 


Substituting  into  Equation  (D-ll) 


N A Sin9 
f r ( © ) h (0 ) 


(D-17) 


This  expression  can  be  evaluated  over  the  image  plane  as  a function  of  0 
to  determine  uniformity  of  focal  plane  illumination  (W) . If  the  pupil  is  elliptical 
the  area  A is: 


A = Trab 


Equation  (D-17)  then  becomes 


Nirab  Sin  9 
fr(9)h(0) 


(D-18) 


(D-19) 
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To  put  this  in  relative  terms,  we  can  refer  to  the  equivalent  F/number  of  the 
nonlinear  lens,  i.e.,  the  circular  aperture  F/number  for  a conventional  lens  that 
yields  the  same  focal  plane  illumination  as  Equation  (D-19). 

If  A was  circular  with  diameter  (D)  and  focal  length  constant  (f)  the  flux 
through  A would  be 


where 


P = NA  — r 

a V 


A = ttD 
4 


(D-20) 


P = NttD2  At 
4S2 


(D-21) 


If  the  lens  is  conventional,  all  this  power  falls  on  an  image  area 

,2 


Therefore , 


« - r 

Ai 


4(fnoe>‘ 


Equating  this  to  Equation  (D-19)  and  solving  for  F. 


ffr(e)h(e) 

abSinG 


(D-22 ) 


This  equation  gives  the  nonlinear  lens  F/number  in  conventional  terms.  It 
is  used  to  assess  F/number  uniformity  of  the  nonlinear  lens  designs. 
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